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Abstract 

We investigate, based on the tight-binding model and in the linear deformation regime, 
the strain dependence of the electronic band structure of phosphorene, exposed to a 
uniaxial strain in one of its principle directions, the normal, the armchair and the zigzag 
directions. We show that the electronic band structure of strained phosphorene, for the 
experimentally accessible carrier densities and the uniaxial strains, is well described by 
a strain-dependent decoupled electron-hole Hamiltonian. Then, employing the decoupled 
Hamiltonian, we consider the strain dependence of the charged-impurity-limited carrier 

mobility in phosphorene, for both types of carriers, arbitrary carrier densities and in both 

* Corresponding author. Tel./fax: -1-98 831 427 4569, Tel: -1-98 831 427 4569. E-mail address: 
yawar. mohammadi@gmail. com 


1 



armchair and zigzag directions. We show that a uniaxial tensile (compressive) strain in 
the normal direction enhances (weakens) the anisotropy of the carrier mobility, while a 
uniaxial strain in the zigzag direction acts inversely. Moreover applying a uniaxial strain in 
the armchair direction is shown to be ineffective on the anisotropy of the carrier mobility. 

These will be explained based on the effects of the strains on the carrier effective masses. 

Keywords: Phosphorene; Tight-binding model; Band Structure; Strain; Carrier Mobility. 

1 Introduction 

Since successful isolation of a single layer of graphite [1] called graphene, as the hrst real two- 
dimensional lattice structure which shows novel appealing properties[2l |3], many researchers 
tried to synthesis or isolate new two-dimensional materials. These efforts resulted in hnding 
other two dimensional materials such as BN[1], transition metal dichalcogenides[5], silicene[6l 
El [HI E] and recently phosphorene. Phosphorene is a single layer of black phosphorus, which 
can be isolated by mechanical exfoliation [101 IH] of black phosphorus. In a single layer of 
black phosphorus, each phosphorus atom covalently couples to three nearest neighbors. This 
conhguration of phosphorus atoms results in a honeycomb-like lattice structure. However, due 
to the sp^ hybridization of s and p atomic orbitals, it forms a puckered surface. The electronic 
band structure of phophorene has been studied using different methods such density functional 
theory calculations [HI EH], k.p method[m [T5] and tight-binding model[T31 ESI EZ]- These 
considerations show phophorene is a direct-band-gap insulator, but with an anisotopic band 
structure. This novel band structure leads to many attractive properties[THl ESI ESI El]- 

Strain tuning is an effective means to tune the physical properties of two dimensional 


2 


materials (for a review, see e.g. Ref. [2l |22]). Puckered structure of phosphorene makes 
this easier, so one can tune and control its electronic and mechanical properties by strain, 
confirmed by recent studies on the effects of uniaxial and biaxial strains in phosphorene 
[l2l|23l[Mll25l[26l|271[2Hll29l[30l[nil3^ These works examined effects of strains applied 
along three principle directions which preserve D 2 h group point symmetry of phosphorene |16] . 
its zigzag and armchair edges and the direction normal to its plane. It has been shown that 
a uniaxial strain in the direction normal can decrease its band gap and even leads to an in¬ 
sulator to metal transition [T21 123 El]- Moreover, the effects of an in-plane uniaxial strains 
along zigzag and armchair edges [281 12^ 130] on the band gap of phosphorene has been stud¬ 
ied. Some other researchers has studied the effects of uniaxial and biaxial strains on the band 
structure [3011311 EH [33] and the optical properties [31] of phosphorene, confirming the capability 
of stain as an effective means to tune the properties of phosphorene. These works showed, when 
the uniaxial strain is applied along the armchair direction, the properties of phophorene change 
further. But, recently, it has been shown [35] that the most effective direction to apply a strain 
and tune the band gap of phosphorene is an in-plane direction (not being along armchair nor 
along zigzag) with a direction angle about 0.2687r counted from the armchair edge. 

According to the high-potential capability of strains to tune the properties of phosphorene, 
driving an analytical relation for the Hamiltonian of strained phosphorene is very desirable, and 
can be used to examine the effects of the strains on the electric, optical and magnetic properties 
of phosphorene. In this paper, starting from the well-known 4-band tight-binging Hamiltonian 
of phophorene [T21ITB] . we obtain a strain-dependent tight-binding Hamiltonian for phosphorene. 
In this paper we work in the linear deformation regime and only consider uniaxial strains. To 
benefit from the D 2 h group point symmetry of phosphorene and reduce the 4-band Hamiltonian 
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to a 2-band Hamiltonian and achieve an analytical result, we restrict our consideration to the 
uniaxial strains applied along three principle directions of phosphorene which preserve its D 2 h 
group point symmetry. Thanks to this symmetry, we can obtain analytical relations for its band 
energies which can be used to explore easily the effects of the uniaxial strains on properties of 
phosphorene. Searching for low-energy structures in strained phosphorene, we use continuum 
approximation and derive the corresponding Hamiltonian dominating low energy excitations. 
Then, by taking into account the weak interband coupling of conduction and valance bands, we 
project the low-energy Hamiltonian into a decoupled Hamiltonian [5BI [STj 155] and show that for 
the experimentally accessible carrier density, the decoupled bands agree well the bands obtained 
from the tight-binding Hamiltonian. Motivated by this fact and the recent studies on the carrier 
mobility in phopsphorene[Tni [20l |33 [30], then we apply our decoupled Hamiltonian to consider 
the strain dependence of the charged-impurity-limited carrier mobility in phosphorene. Our 
result shows that one can tune the amount and the anisotropy of the mobility in phosphorene, 
by making use of a uniaxial strain in the normal or zigzag direction. 

The rest of this paper is organized as follows. In Sec. H we reproduce the known 4- 
band Hamiltonian of phosphorene. In Sec. HI we explain how one can, in general, insert the 
effects of the strains in the Hamiltonian and obtain a general formalism for the strain-dependent 
Hamiltonian of phosphorene. Sec. IV devoted to consider the strain dependence of the charged- 
impurity-limited carrier mobility in phosphorene. We end the paper by summarizing our results 
in Sec. V. 
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2 Structure and tight-binding Hamiltonian of phospho- 


rene 

The lattice structure of phosphorene and the necessary lattice parameters to construct the 
tight-binding Hamiltonian of phosphorene, including the lattice constant, the bond angles and 
the transfer energies, have been introduced in Fig. [H The unit cell of phosphorene (solid-line 
rectangle in Fig. [1]) consists of four phophorus atoms, two atoms in the lower layer represented 
by the grey circles (called A and B) and two atoms in upper layer represented by the red circles 
(called C and D). Hence, the tight-binding Hamiltonian of phophorene can be written in terms 
of a 4 X 4 matrix as 


0 

tylB(k) 

^^c'(k) 

tylD(k) 

^BA(k) 

0 

^Bc(k) 

^BU>(k) 

tcA^X) 

^CB(k) 

0 

^CD(k) 

^DA(k) 


^Dc(k) 

0 


acting in (0^, -ipB, i’Ci '^PdY' with k being the two-dimensional momentum. Notice that = 

Moreover, it has been shown [T31 [TB] that if we only retain the transfer energies up to 
the hfth nearest neighbors, the tight-binding approximated band structure agrees well with 
its density functional theory band structure. These transfer energies are|l3] ti = —1.220 eV, 
t 2 = -1-3.665 eH, ts = —0.205 eH, = —0.105 eV and = —0.055 eV. So we can rewrite the 
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Hamiltonian matrix as 


H. = 


/ik + /: 


3k 


0 


f* f* f* 

i4k i2k “T Jbk 

f* f* f* 

i2k ^ iSk Ji 


Uk 


/ 4 k 

/ 2 k + / 5 k 

/ 2 k + / 5 k 

/ 4 k 

0 

flk + fsk 

flk + fsk 

0 

cos(fcy|/i), / 2 k = 

Here 

h) = {xi, 


( 2 ) 


is drawn from A (The origin of the cartesian coordinate system) to one of the ith near¬ 
est neighbors (See Fig. [T]). They are fi = (—di cos a, di sin a, 0), r 2 = (d 2 cosd, 0, d 2 sin d), 
fh = (di cosa-f 2 d 2 cosd, di sina, 0), = (—di cos a — d 2 cos d, di sina, d 2 sin d) and rs = 

(—2diCosa — d 2 cosd, 0, d 2 sind) where cosd = —One can take into account the D 2 h 
group point symmetry in phosphorene and project the four-band Hamiltonian into a reduced 
two-band Hamiltonian as[ 

/ 


= 


\ 

/ 


(3) 


/ 4 k flk + / 2 k + /sk + /sk 

flk + / 2 k + fsk + fbk hk 

acting in (0^ -|- (pc, <Pb + (po )^The corresponding energy bands, obtained by diagonalizing 
the Hamiltonian matrix, are given by 


Ek — / 4 k ± l/lk + / 2 k + /sk + /skl, (4) 

where +(—) denotes to the conduction(valance) band. We have shown the energy spectrum 
of phosphorene obtained from two-band Hamiltonian in Fig. [2j It is evident that minimum 
(maximum) of the conduction (valance) energy band is at F point. If we apply continuum 
approximation to the obtained two-band Hamiltonian and retain the terms up to the second 
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order in k, we can reproduce the known Hamiltonian of phosphorene |1 6 [ [T7] . 


Hu = 


u + r]^kl + r]ykl 5 + '^^kl + '^yk^ + ixkx ^ 


( 5 ) 


^ ^ + Ixkl + Xykl - ixK U + r]xkl + r^ykl ^ 

where u = At^ = 0.42 eH, r]x = — 2 ^ 4 X 4 = 1.03 rjy = —2t/^y1 = 0.56 5 = 2ti + ^2 + 

2^3 + ^5 = 0.76 eH, 7 ^; = —tixf — — tsxl — yXg = 3.51 el/A^, 7 ^^ = —hyl — = 3.81 eV 

and X = 2tiXi + 12 X 2 + 2 ^ 3 X 3 + ^ 5 X 5 = —5.34 cHA which agree well with the other calculations 
[I?] (Notice that in our calculations the zigzag edge lies along the x-axis.). The corresponding 
energy spectrums are given by 


Eu = u + Vxkl + %kl ± ^((5 + Xxkl + XyklY + 


( 6 ) 


where +(—) denotes to the conduction (valance) band. It is evident that the energy spectrum 
is linear in the kx direction while in the ky direction it is parabolic. Due to the large band gap, 
which leads to a weak interbans coupling, one can decouple the electron and the hole bands 
into a low energy regime. In this approximation, Eq. [6] can be written as[r71 [ 3 S 1 EZl ES] 


21,2 


u + r]xkl, + rjyky ± 5(1 + 2 


( 7 ) 


In this approximation the electron and hole effective masses in the x and y directions are given 

^ex = 2(,.+J+xV 2^) = 0-168 mo, mey = = 0.852 mo, mux = 2 ( 7 .-J-xV 2 ^) = 

0.184 mo and rrihy = 2 ('y^-y ) “ 1-146 mo, which mo is the mass of a free electron, in good 
agreement with recent result [T3] . To see that in what region this approximated energy bands 
agree well with the other results, we have shown all three set energy bands obtained from the 
tight-binding, the low energy and the decoupled Hamiltonians, in Fig. [2j One can see that 


in the ky direction all three set energy bands agree well in a wide range of the energy and 
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the momentum. Moreover, in the kx direction the tight-binding energy bands and the low- 
energy bands agree well too, but the decoupled bands overlap with them only up to 0.14 eV 
(0.13 eV) with respect to the bottom (top) of the conduction (valance) bands. This corresponds 
to n = 2.20 X 10^^ and n = 2.44 x 10^^ for the electron and hole densities. These 
indicate that the low energy excitations in phosphorene are well described by the decoupled 
Hamiltonian [13 IMl EU 138] . 

3 Strain-dependent tight-binding Hamiltonian 

In this section we rederive the tight-binding Hamiltonian of phosphorene in the presence of the 
uniaxial strains applied along the principle directions of phosphorene. To insert the effects of 
the applied strain in the tight-binding Hamiltonian of phosphorene, hrst we must determine the 
effects of the strain on the transfer energies and the bond lengths. It has been shown [39] that the 
transfer energies between s and p orbitals, which construct the electronic bands of phosphorene, 
depend on the bond length as t oc To obtain this relation, it has been supposed that the 
applied strain doesn’t change the bond angles and only change the bond lengths. Within the 
linear deformation regime, this is a reasonable assumption. Since the change in the bond angles 
in a strained lattice, at least, includes the terms of second order in terms of the applied strain. 
So they can be ignored in the linear deformation regime. Hence, we only need to determine the 
strain dependence of the bond lengths and the other inter-atomic distances. 

Let us construct our formalism in a general case in which phosphorene is exposed to strains 
applied along all three principle directions of phosphorene, the armchair (x-direction) and the 
zigzag (y-direction) edges and the the direction normal to the phosphorene plane (z-direction). 


So, the deformed coordinates are given by 




X 


y 

V ^ / 


/ 


\ / \ 

X 


V 




y 

V ^ / 


( 8 ) 


1 + ea; 0 0 

0 1 + ey 0 

0 0 1 + 

where ex, ey and are the normal strains applied along the x-, y- and z-directions respectively. 
In this paper, we restrict onr considerations to the linear deformation regime, so the bond 
lengths and the other atomic distances, in general, can be expanded in terms of all components. 


ex, ey and ez as 


= r + axex + Uyey + 


(9) 


where r and = \J[x'^Y + f^e unstrained and strained bond lengths respec¬ 

tively and ttx, Ciy and are the strain-related geometrical coefficients, given by ax = ^|erj=o = 
ay = ^ and If we insert Eq. |9]into the relation t cc and 

expand it in terms of the strains and only retain the terms up to the hrst order in e we get 


f = t --{axex + ayey + azez)t, (10) 

where t and are the unstrained and the strained transfer energies respectively. 

As mentioned above, the strains applied along all three principle directions of phosphorene 
don’t break D 2 h symmetry. So the electronic excitations in strained phosphorene are dominated 
by the reduced two-band Hamiltonian, Eq. [3l but after substituting the deformed transfer 
energies and bond lengths into it. Recalling the relations obtained for the strained bond lengths 
and the transfer energies, and substituting them into the tow-band Hamiltonian, Eq. [3l we get 


Hv = 


L 


4k 


fit + /Ik + /Ik + fl 


/ik + /Ik + /Ik + /Ik ^ 


e* 

5k 


/ 


4k 


( 11 ) 


9 











where /ik = cos{kyyl), / 2 k = / 3 k = cos{kyyl), / 4 k = cos(/c^a;^) cos{kyyl) 

and /sk = The corresponding electron and hole energy bands are given by 


-^k — /4k ± l/ik + /2k + /4 + /skl- 


( 12 ) 


with +(—) denoting to the electron (hole) band. It easy to show that strained phopsphorene 
has a direct band gap at T point, in agreement with the recent density functional theory 
[T 2 II 23 [3ll|28] and tight-binding|10] calculations done in the linear deformation regime. Similar 
to the unstrained case, to capture the low energy physics of strained phosphorene, one can 
expand the matrix elements around T point and only retain the terms up to second order in k 
and first order in e. Hence, the low energy Hamiltonian becomes 


Hk = 


^ u‘ + 7]lkl + Tilkl 6'- + YM + + itK '' 


(13) 


lx X ' ly y ' lx X ' ly y 

- ix^k^ + 7]lkl + ylkl j 

The values of the matrix elements depends on the directions in which the strains are applied 

(See appendix). The applied strains can affect, by changing the energy gap and on the 

interband coupling. When the interband coupling is weak, one can project the low energy 

two-band Hamiltonian into a decoupled Hamiltonian which is given by 

\ 


Hk = 


2m% 


2m% 


\ 


0 


0 




where El = u’' + 5^^, E\ = — S'" and 


2ml 




(14) 




2 (hx + 7S + (x^)V 2 ^" 


= TT, -T 




= 


^hy = 


2(7x -Vx + (x")V 2 ^ 0 ) 


2 ( 7 ^-7^)' 


(15) 
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In the reminder of this section we consider the effects of the strains applied along all three 
principle directions of phosphorene on its electronic band structure, as a key feature of crys¬ 
talline materials to explore their other physical properties. We have two aims. One is to 
see whether our tight-binding Hamiltonian reproduces previous results [T^ [271 IMl l28l HO] for 
the energy gap of strained phosphorene. The other is to show that in what energy region 
the decoupled energy bands agree well with the others, obtained from the low-energy and the 
tight-binding Hamiltonian. 

Uniaxial strain along the normal direction (z-axis)- Let us hrst explore effects of a uniaxial 
strain in the normal direction (z-direction), e^, = = 0 and 7 ^ 0. If we recall the relations 

obtained for the the strained bond lengths and the transfer energies, and substitute them into 
Eq. [T 2 I we get 


AE, = 


— ( 8 fi—-|- 4 f 2 ~|' + 8 ^ 3 —-|- 4 ^ 5 —|■)e^ — —12.693e^, 


(16) 


for the strain-induced modulation in the energy gap. This shows that the energy gap of phos- 
phorenre decreases (increases) linearly when it is exposed to a uniaxial tensile (compressive) 
strain in the normal direction. This is in agreement with the previous hrst-principle [T^ l27l l34] 
and tight-binding | 1 D] studies on the strain-induced modulation in the energy gap of strained 
phosphorene, done in the linear deformation regime. This can also be seen in Fig. [3] where we 
have shown the energy bands of strained phosphorene for different values of obtained by di¬ 
agonalizing Eq. [TT] (black curves), Eq. [12] (red dashed curves) and Eq. [TH (green dotted-dashed 
curves). In this hgure right (left) panels show the energy bands of phosphorene in the presence 
of a uniaxial tensile (compressive) strain applied in the normal direction, and in each panel the 
energy bands have been drown in both T — X and T — E directions. 
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This figure also shows that a uniaxial tensile (compressive) strain in the normal direction 
enhances (weakens) the anisotropy of the both electron and hole energy bands slightly (See 
black cures in Fig. [3]). This becomes more clear in the next section, where we consider the 
strain dependence of the carrier mobility in strained phosphorene. Moreover one can see that in 
the presence of a uniaxial tensile (compressive) strain, the energy range in which the decoupled 
bands agree well with the other bands becomes limited (extended). This is mainly due to the 
effect of the stain on the energy gap (See Eq. [16]). When the energy gap increases (decreases), 
the coupling of the conduction and the valance band is enhanced (weakened) and the decoupling- 
band approximation becomes more (less) accurate. For a uniaxial tensile (compressive) strain 
about ez = 0.04 (e^ = —0.04), the decoupled conduction band overlap well with the tight- 
banding conduction band up to 0.11 eV (0.18 eV) with respect to the bottom of the conduction 
band. By making use of n = where Ep is counted from the bottom (top) of the 

conduction (valance) band, one can show this agreement corresponds to n = 1.30 x 10^^ 

(n = 3.25 X 10^^ cm“^) electron density. This agreement for the valance band is up to 0.10 eV 
(0.16 eV) with respect to the top of the valance band, corresponding to n = 1.48 x 10^^ cm~^ 
{n = 3.60 X 10^^ cm~‘^) hole density. 

Uniaxial strain along the zigzag edge (y-axis)- When phosphorene is exposed to a uniaxial 
strain along its zigzag edge, the strain-induced modulation in its energy gap is given by 

AEg = — (8fi— 4^2—I + 8^3— I-|-4^5— |)ej, = 5.945ej^, (17) 

ri r2 r3 

which shows that a uniaxial tensile (compressive) strain in the zigzag edge increases (decreases) 
linearly the energy gap. This agrees well with the recent studies [T0ll28ll40] . Figure Hjshows that 
in the presence of a uniaxial tensile (compressive) strain along the zigzag edge, the anisotropy 
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of the band structure is weakened (enhanced). Moreover it is evident that, in the presence of 
a uniaxial tensile (compressive) strain about ey = 0.04 {ey = —0.04), there is good agreement 
between the decoupled and the tight-binding conduction bands up to 0.16 eV (0.13 eV) with 
respect to the bottom of the conduction band. In the valance band the overlapping is up to 
0.15 eV (0.12 eV) with respect to the top of the valance band. 

Uniaxial strain along the armchair edge (x-axis)- In the presence of a uniaxial along the 
armchair edge of phosphorene, the strain-induced modulation in its energy gap is given by 


^Ey = 


-(8fi4 + 4^2^ + 8^3^ + 4t5^)ex = 3.708e, 

7^2 ^3 "^5 


(18) 


which shows that the energy gap is a linear function of the applied strain, and increases (de¬ 
creases) when phosphorene is exposed to a uniaxial tensile (compressive) strain in agreement 
with the recent studies [T0| 140] . Comparison of Eqs. ITTI and flSl shows that, for same uniaxial 
strains along the zigzag and armchair edges, the uniaxial strain along the zigzag edge induces a 
larger band gap variation. Effects of the applied strain on the anisotropy of the band structure 
can be seen in Fig. |0] which, as it is expected, is unlike the effects of the uniaxial strain along 
the armchair edges. In the presence of a uniaxial tensile strain about Cy = 0.04 along the zigzag 
edge, the overlapping of the decoupled band with the tight-binding is up 0.16 eV and 0.14 eV 
for the conduction and the valance bands respectively, while for a compressive strain about 
Cy = —0.04 they agree only up 0.13 eV for both conduction and valance bands. 

We end this section by this conclusion that the electronic band structure of strained phos¬ 
phorene, for the experimentally accessible carrier densities and the uniaxial strains applied along 
all three principle directions of phosphorene, is well described by the decoupled Hamiltonian. 
Motivated by this fact, we apply it to consider strain engineering the charged-impurity-limited 
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carrier mobility in phosphorene. 


4 Strain engineering the charged-impurity-limited car¬ 
rier mobility 


In this section, employing our strain-dependent decoupled Hamiltonian, we investigate the 
strain dependence of the impurity-limited carrier mobility in phosphorene for both types of 
carriers, electron and hole, and along both armchair and zigzag edges. The carrier mobility, 
/i, is dehned as fx = ajne where a is the electrical conductivity, n is the carrier density and 
e is the electron charge. To calculate the electrical conductivity we use the semi-classical 
Boltzmann transport theory combined with the relaxation time approximation. Moreover we 
restrict our calculation to the steady state and suppose that the two-dimensional electron gas 
in phosphorene is homogenous, so the electrical conductivity is given by 



where i is x, y, Qs = 2 is the spin degeneracy, k = (k^, ky) is the two-dimensional momentum and 
Vi = hki/mj is the electron velocity in the i direction with ki and mj being the corresponding 
electron or hole momentum and mass. i?k is the energy band obtained from the strained- 
dependent decoupled Hamiltonian (Notice we have omitted the electron and hole indexes in 
ml and E]^.), /(i?k) is the Fermi-Dirac distribution function and T(i?k) is the relaxation time. 
Let us suppose that the impurities are static, of symmetric potential and have no internal 
excitations. So the relaxation time is given by 


1 27rni 

t{e^) ^ ~ir 


{‘2t^Y e(g) 


2(l-cosV)d(^fc 



( 20 ) 
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where rij is the number of impurities per unit area, q = |k — k'| and is the scattering 
angle between k and k\ Vi{q) = is the Fourier transform of the potential of the charge 
impurity and k = {Ksub + i^enc)/‘^ is the effective dielectric constant with Kgub and Kenc being the 
dielectric constant of the substrate {i^sub = 2.5 for Si02 substrate [19]) and the encapsulating 
layer respectively which for vacuum is zero. e{q) is the dielectric function which within the 
random phase approximation is given by e{q) = 1 + ^^n(g), where n(g) is the polarizability 
function. The polarizability function can be written [19] as 


n(g) = 


mi.m\ 


X -y 


Trr 




( 21 ) 


^ (ll/ml + Qy/^l 

where Ep is the Fermi energy of strained phosphorene for a fixed carrier concentration. 

If we introduce new variables asp^ = andpy = we have Ep = —(Px+ 

X “'y '‘'eff 

Py) for the energy bands with mljj = ^m%my. In the new momentum, space the electrical 
conductivity is given by 

df{Ep) 


— 


mf I 


dEr, 




( 22 ) 


leading to an = for the electrical conductivity of strained phosphorene at 


zero 


(p^ 

temperature, where Ep = „ f with pp being Fermi momentum in the new momentum space. 


e// 


t{Ep) is given by 


riim 


eff 


r(Ei.) 


TT?r Jo 


dd 


27re^ 


y/^KkpVJ — COS 6 + 2ne‘^D{Ep) 


(1 — cos6'). 


(23) 


where D{Ep) = is the carrier density of states at the Fermi energy and 


kp — 


\ 


27m 


TTl^ 771^ 

(^)i/2 cos2 0+ (^)V2sin2 0 

m% ml 


(24) 


is the anisotropic Fermi momentum. In Eq. [2ll 9 is counted from the x-axis, and n = JJ^Ep is 
the carrier density in strained phosphorene, being a linear function of Fermi energy as same as 
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the carrier density in the ordinary two-dimensional electron gas. Hence the zero-temperature 


carrier mobility in strained phosphorene is given hjfiu = 


eriEi,) 


In Fig. [6] we have shown our numerical results for the strain dependence of the charged- 
impurity-limited electron (left panels) and hole (right panels) in phosphorene exposed to the 
uniaxial strains in the normal direction (z-axis). The upper (lower) panels shows the carrier 
mobility along its armchair (zigzag) edge, and orange {n = 0.2 x 10^^cm“^) to black [n = 
1.0 X 10^^cm“^) curves correspond to different carrier densities with An = 0.2 x lO^^cm”^. The 
density of the charged impurities is supposed to be rij = 1.0 x 10^^cm“^ which is typical of the 
Si02 substrate. Figure [6] shows that the carrier mobility along the armchair direction is higher 
than that along the zigzag direction, as same as that in unstrained phosphorene [TOl 1^ 13711^ . 
This can be understood by this fact that, in the presence of both uniaxial tensile and compressive 
strains, the carrier effective mass along the armchair edge is always smaller than that along 
the zigzag edge. This can be tested by making use of the Eqs. [151 [251 and [251 Moreover 
one can see that the carrier mobility along both armchair and zigzag directions increases by 
increasing the carrier density. This is the familiar feature of the ordinary two-dimensional 
electron gas^T], arising from the linear dependence of its carrier density on the Fermi energy 
(In phosphorene the carrier density depends on the Fermi energy as n = Figure [6l 

also shows that in the presence of a tensile (compressive) strain in the normal direction, the 
carrier mobility along the armchair edge increases (decreases), while the carrier mobility along 
the zigzag edge decreases (increases). This property originates from the effect of the strain on 
the anisotropy (and consequently the carrier effective mass) in phosphorene, as explained in the 
previous section. To explain this property further, we rewrite the relation of the carrier mobility 


flu — 6 


it rnl^^ 


. It is easy to show that the effect of the strain on part is weak and 

e// 
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it mainly affects on part. By making use of the Eqs. [TSl [23 and[2Hl one can show that 

‘‘'i 

applying a uniaxial tensile (compressive) strain in the normal direction decreases (increases) 
both electron and hole effective masses in the armchair (zigzag) direction, and consequently 
their mobilities in the armchair (zigzag) direction increase (decrease). 

In Fig. [7|we have compared the effect of the direction of the applied straina on the carrier 
mobility in phosphorene. Figure [7| shows that, unlike the strains in the normal direction, 
applying a uniaxial tensile (compressive) strain in the zigzag direction decreases (increases) 
both electron and hole mobilities in the armchair (zigzag) direction. This originates from their 
different effects on the anisotropy (and consequently the carrier effective mass) in phosphorene, 
as explained in the previous section and above. This hgure also shows that applying a uniaxial 
strain in the armchair direction weakly change the carrier mobility in phosphorene. Moreover 
Figs. [6] and [7| show that applying a uniaxial tensile (compressive) strain in the normal (zigzag) 
direction enhances (weakens) the anisotropy of the carrier mobility in phosphorene. While in 
the presence of a uniaxial compressive (tensile) strain in the normal (zigzag) direction, the 
carrier mobility is weakened (enhanced). 

5 Summary and conclusions 

In Summary, we investigated the electronic band structure of strained phosphorene within the 
linear deformation regime and based on the tight-binding model. We restricted our consider¬ 
ation to the uniaxial strains applied along one of the principle directions of phosphorene, the 
normal, the armchair and the zigzag directions. We showed that the derived strain-dependent 
energy spectrums reproduce the previous results for the energy gap of strained phosphorene. 
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Then we applied the continuum approximation to derive the corresponding low-energy Hamilto¬ 
nian. Moreover we showed, when the interband coupling is weak, the low-energy Hamiltonian 
can project into a decoupled electron-hole Hamiltonian. We found that the electronic band 
structure of the strained phosphorene, for the experimentally accessible carrier densities and 
the mechanical strains, is well described by the decoupled Hamiltonian. Motivated by this fact 
we used our strain-dependent decoupled Hamiltonian to investigate the strain dependence of 
the charged-impurity-limited carrier mobility in phosphorene. We examined the dependence of 
carrier mobility on the direction of mobility, the carrier type, the carrier density and the direc¬ 
tion of the applied strain. We showed the dependence of the carrier mobility on the direction 
of mobility, the carrier type and the carrier density is same as that in unstrained phosphorene. 
Moreover, as a point worthy of mention, we found that applying a uniaxial tensile (compressive) 
strain in the normal direction decreases (increases) carrier mobility in the armchair (zigzag) 
direction. While in the presence of a uniaxial tensile (compressive) strain in the zigzag direc¬ 
tion the carrier mobility is decreased (increased). We also showed that a uniaxial strain in the 
armchair direction don’t changed the carrier mobility approximately. These properties were 
explained based on the effect of the applied strain on the anisotropy of the carrier effective 
mass in phosphorene. 
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A Appendix: Calculating the elements of the strain- 


dependent low-energy Hamiltonian matrix 


The matrix element in Eq. [131 the linear deformation regime, depend in general on the 
applied strain as 

u = u + eu 

vl = Vx + ehx 

nl = ny + (^Vy 
= 5 + eS' 

7^ = 7x + eql 
7^ = 7y + n'y 

x" = X + ex, (25) 


where the coefficients of the applied strain, for a uniaxial strain in the normal direction, are 
given by 


u = 


Vx 

% 

6 ' 

lx 

% 

/ 

X 


= 


-8^44 = 0.32 eE 
rf 

^2 2 

4^4 = —0.75 eV 
ri 

4t4^ = -0.43 eVA^ 
ri 

-4ti^ - 2t2% — - 2t5% = -6.58 eV 

Ti 


^2 


^3 


r‘2 

^5 


= 2ti- 


^2 2 

r'i 

'1 

v 2„,2 


+ t 


^2 2 
^2-^2 

'2 


+ 2t; 


^2 2 


+ h 


^2 2 

= 1.44 


2„,2 


2ti^ + 2^3^ = 0.00 el/A^ 


= -4ti 


zfxi 

r‘2 

'1 


- 2tl-k^ = 4.74.1/A. 


r'2 

'2 


ri 


(26) 
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These coefficients, when the strain is applied in the zigzag edge (y-axis), become 


u 


V 


X 


% 

6 ' 

lx 


1 


y 


/ 

X 


-8t4^ 

= 0.19 el/ 


rj 



it/A 

= -0.45 eVA^ 


rj 



AtJ^- 

4^41/4 = 0.89 eVA^ 

ri 




-2t/4-4t/4 

2t ^5 

rf 

rf rf 

r| 

yWi 

III 2 

VnXn Vixi 1 

+ t2 „ + Zt3 0+15 

rf 

^2 

r-s 

2«i4- 

2t,yf + 2h4- 

2hyl = 1 

rf 

ri 


ry^ ry^ 


.2^2 


^2 

' 5 


= -3.81 eVA^ 


2 t 5 ^ = 3.43 eVA, 
ri 


(27) 


while for a uniaxial strain along the armchair edge, they are 


u 


V 


X 


% 

S' 


1 


X 


1 


y 


/ 

X 


^2 

-8^4^ 

ri 


= 0.33 el/ 


At 


4^4 


XI 

xlvl 


Atixl = 1.21 eV A^ 
= -0.45 el/A^ 


X 


X 


Xn 


Xt 


-AU - 2to^ - AU^ - = 2.15 el/ 


2 

rf 

xf 

2ti^ - 
rf 

xfyl 


r2 


4 

X 2 


n n 

2 


2tiX-^ + ~ ^23^2 ffi 2^3 


ylA 


r2 
' 2 


2^3X3 + ^5 


yi4 


hxl 


2ti 


r 


+ 2^3 


1 

xf 

-4ti^ 

rf 


4yi 

rj 


= -3.81 el/A^ 


+ 2tiXi — 2t 


2 2 

ri 


+ t 2 X 2 — 4^3— + 2^3X3 — 2^5— + ^5X5 

ri 


r2 

^5 


2.37 el/A^ 


-2.96 eVA. (28) 


The lattice parameters, which we used here, are di = 2.22A, d 2 = 2.24A, ti = —1.220eV, 
t 2 = 3.665el/, ts = -0.205el/, U = -0.105el/, h = -0.055el/ [I3], a = 0.26757r and 9 = 0.5677r 

[3S]. 
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Figure 1: (a) The top view of phosphorene lattice structure, ti indicates to transfer energies 
from a site to its ith nearest neighbors. The solid-line (dashed-line) rectangle denotes to the 
unit cell in the 4-band ( 2 -band) model. The other parameters are di = 2.22A, ^2 = 2.24A,[T3] 
and a = 0.26757r.[35] (b) The side view of phosphorene lattice structure where (3 = 0.5677r.|35] 
(c) The coordinate system used in this work. The armchair edge is supposed to be along the 
x-axis and the zigzag edge along the y-axis. 
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Figure 2: The energy spectrums of phosphorenere around F point obtained from the two-band 
tight binding Hamiltonian (solid black cure), the low-energy Hamiltonian (red dashed curve) 
and the decoupled Hamiltonian (green dotted-dashed curve). 
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Figure 3: The energy spectrums of strained phosphorenere around F point obtained from the 
two-band tight binding Hamiltonian (black solid cure), the low-energy Hamiltonian (red dashed 
curve) and the decoupled electron-hole Hamiltonian (green dotted-dashed curve). Phosphorene 
is exposed to uniaxial tensile (left panels) and compressive (right panels) strains applied in the 
normal direction. 
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Figure 4: Same as Fig. |3]but for strains applied along the zigzag edge. 
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Figure 5: Same as Fig. |3]but for strains applied along the armchair edge. 
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Figure 6: The strain dependence of the charged-impurity-limited electron (left panels) and 
hole (right panels) mobilities along the armchair (upper panels) and zigzag (lower panels) 
edges of phosphorene exposed to uniaxial strains in the normal direction. Orange to black 
lines correspond to n = 0.2 x to n = 1.0 x carrier densities with An = 

0.2 X and the density of the charged impurities is rij = 1.0 x lO^^cm”^. 


30 










- 3 - 2-1 0 1 2 3 4 

strain(%) 


Figure 7: The strain dependence of the charged-impurity-limited electron (left panels) and hole 
(right panels) mobilities along the armchair(upper panels) and zigzag(lower panels) edges of 
phosphorene, when it is exposed to uniaxial strains along the normal (black cure), armchair 
(red cure) and zigzag (green cure) directions. The electron and hole densities are 1.0 x 
and the density of the charged impurities is n* = 1.0 x 10^^cm“^. 
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